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ABSTRACT 

We investigate the signature of primordial non-Gaussianities in the weak lensing bispectrum, 
in particular the signals generated by local, orthogonal and equilateral non-Gaussianities. The 
questions we address include the signal-to-noise ratio generated in the Euclid weak lens- 
ing survey (we find the Icr-errors for /ml are 200, 575 and 1628 for local, orthogonal and 
equilateral non-Gaussianities, respectively), misestimations of /nl if one chooses the wrong 
non-Gaussianity model (misestimations by up to a factor of ±3 in /nl are possible, depend- 
ing on the choice of the model), the probability of noticing such a mistake (improbably large 
values for the ;t'^ -functional occur from /nl ~ 200 on), degeneracies of the primordial bispec- 
trum with other cosmological parameters (only the matter density Q.„, plays a significant role), 
and the subtraction of the much larger, structure-formation generated bispectrum. If a prior 
on a standard wCDM-parameter set is available from Euclid and Planck, the structure forma- 
tion bispectrum can be predicted accurately enough for subtraction, and any residual structure 
formation bispectrum would influence the estimation of /nl to a minor degree. Configuration- 
space integrations which appear in the evaluation of ;^'^-functionals and related quantities can 
be carried out very efficiently with Monte-Carlo techniques, which reduce the complexity by a 
factor of (9(10'*) while delivering sub-percent accuracies. Weak lensing probes smaller scales 
than the CMB and hence provide an additional constraint on non-Gaussianities, even though 
they are not as sensitive to primordial non-Gaussianities as the CMB. 
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1 INTRODUCTION 

As cosmological data improves, it is becoming increasingly feasi- 
ble to probe models of the early universe. In particular, primordial 
non-Gaussianity has emerged as a leading window on to the physics 
of inflation and the early universe ( for reviews, see iBartolo et al.l 
l2004lKomatsull2010l : lLangloisll201ll) . Although non-Gaussian per- 
turbations could in principle take any form, in practice searching 
for just a few shapes of the bispectrum (3-point correlation func- 
tio n) allows one to discriminate between entire classes of mod- 
els I Bernardeau & Uzanl2003l : lKomatsu et alj2009l : lKomatsu et al.l 



and gravitational lensing I Fedeli et alj l201ll ; IPace et al.l 1201 ll : 
iMarian et alj201 lUjeong et alj201 ih 

Most constraints on non-Gaussianities are reported using 
CMB-observations, either by m easuring the bispectrum of the 



UOll h. Many alternative models to inflation produce the same non- 
Gaussian shapes, and therefore can also be constrained at the same 
time. 

The theory of primordial non-Gaussianities is very evolved, 
and predictions for non-Gaussianities in different observational 
channels have been made: number counts and large-scale structure 
statis tics (Verde 2010; Desjacques & Seljak 2010a, b; Fedeli et al. 
l20Ilh . the cos mic microwave background even outside the Sachs- 
Wolfe regime l IFergusson & ShellarduOOTI : iFergusson et al.ll2010f) 
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temperature perturbation directly (Komatsu 2003; Creminelli et al 
2006 ; Yadav & Wand ell 2008; Curt oetal. 2009; Ko matsu etal 
20091 ; IVielva & Sand |2009| . l20ld ; ICasaponsa et al.l 1201 ll) , by 
measuring the skewness of weighted averaged CMB -patches 
JMukheriee & Wand |2004) or by quantifying the correspond- 
ing Minkowski functional JCabella etal] 120051 ; [Oott et alj|2007l ; 
JHikage et al. 2008). The tightest bounds on the ampl itude /nl of the 
bispec trum, -10 < /nl < 74 has been obtained bv lKomatsu et al.l 
( l201lh . 

So far, observational constraints on the weak lensing bispec- 
trum JBemardeau et al.l2003l ; IShi et al.l201 1 ) mainly concerned the 
non-Gaussianities generated b y structure formatio n and their break- 
ing of the n„,-o"8 degeneracy JMenard et alj|2003l) . and first results 
have been obtained using the skewness of the aperture-mass statis- 
tic ISchneider et al. 1998; Semboloni et al. 2011). 

In this paper we forecast the constraints which the weak lens- 
ing bispectrum will be able to make on primordial non-Gaussianity, 
especially with a view to Euclid. Although the constraints are not 
competitive with the CMB if the non-Gaussianities are scale inde- 
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pendent, they are complementary since they probe smaller scales 
compared to the CMB and provide constraints on a possible scale- 
dependence I Lo Verde et alj|2008h . A scale dependence of /nl at 
the same order as the s pectral index of the power spectrum is natu- 
ral JBvmes et alJuOlOf) and it may be much stronger, e.g. a "step- 
function" which is ze ro on large scales and large on small scales 
JRiotto & SlothluOllI) . Weak lensing also has lower systematic er- 
rors than other large scale structure probes such as the galaxy bis- 
pectrum, scale dependent bias and cluster counts. Because the weak 
shear provides a linear mapping of the cosmic matter distribution, 
the statistical properties of the source field are conserved in the ob- 
servable. 

After a brief recapitulation of cosmology and structure forma- 
tion in Sect.|2]we introduce primordial and structure formation non- 
Gaussianities in Sect. [3] in particular the bispectral shapes and the 
motivation for studying them. The mapping of non-Gaussianities 
by weak gravitational lensing is treated in Sect.|4] where we inves- 
tigate the properties of the weak lensing bispectrum in its scale- and 
configuration dependence, and how it builds up as a function of sur- 
vey depth. Statistical questions concerning the signal strength and 
misestimations of the non-Gaussianity parameter are addressed in 
Sect.[5] before we focus on systematic errors in the non-Gaussianity 
parameter due to incompletely removed structure formation non- 
Gaussianities in Sect.|6] We summarise our main results in Sect. [7] 
and provide visualisations of the weak lensing bispectrum sourced 
by different non-Gaussianity shapes in Appendix IaI 

The reference cosmological model used is a spatially flat 
wCDM cosmology with Gaussian adiabatic initial perturbations for 
the cold dark matter density. The specific parameter choices are 
n„, = 0.25, tis = 1,0-8= 0.8, Q.,, = 0.04 and Hq = 100/ikm/s/Mpc, 
with h = 0.72. The dark energy equation of state is constant in time 
with a value of w = -0.9. 



2 COSMOLOGY AND STRUCTURE FORMATION 

In spatially flat dark energy cosmologies with the matter density 
parameter Cl,„, the Hubble function H(a) = dlna/dt is given by 



HHa) 



//2 ~ ~a^'^ a3(l+vr) ' 



(1) 



for a constant dark energy equation of state-parameter w. Comov- 
ing distance x and scale factor a are related by 



X- 



r' da 

■j„ a^H{a)' 



(2) 



such that the comoving distance is given in units of the Hubble 
distance xh = c/Hg. For the linear matter power spectrum P(k) 
which describes the Gaussian fluctuation properties of the linearly 
evolving density field S, 



(S(k)S(k')} = {IjifSoik + k')P(k) 



(3) 



the ansatz P(k) oc k"' T-(k) is chosen with the transfer function T(k), 
which is well approximated by the fitting formula 



(4) 



ln(1.2.34,) 
^ 2.34^ '^^ 

for low-matter density cosmologies l IBardeen et al.l 119860 . The 
polynomial /7(^) is given by p(^) = l + 3.89^-l-(16.1^)^-l-(5.46(?)''-l- 
(6.71^)'*. The wave vector k = qT enters rescaled by the shape pa- 
rameter r JSugivamall 1 9951) . 



r = Q.i„h exp 



-Q,,\l + 



V2A Y 

a, j 



(5) 



The fluctuation amplitude is normalised to the value ag on the scale 
R = 8 Mpc/h, 



:^ jdkk'wl(k)P(k), 



(6) 



with a Fourier-transformed spherical top-hat W^ik) = 3ji(kR)/(kR) 
as the filter function. je(x) denotes th e spherical Bessel function 
of the first kind of order £ l lAbramowitz & Stegun 1972). The lin- 
ear growth of the density field, d(x,a) = D+(a)S(x,a = 1), is de- 
scribed by the growt h function D+(a), which is the solution to the 
growth equation (Tu rner & Whitelll997l : lwang & Steinhardtll 19981: 
JLinder & Jenkinsll2od3[K 

d^ 1 / dlnH\ d 3 

-^D+(fl) + -\3 + — —£•+(«) = -^Q„,{a)DAa). (7) 

da-^ a \ dlna / da 2a^ 



3 NON-GAUSSIANITIES 

3.1 Primordial non-Gaussianities 

We write the primordial bispectra in te rms of the Ba rdeen curvature 
perturbation O cBardeen.l980. ; .Bardeen et al.lll983l) , which may be 
related to th e primordial curvature perturbation f = SO/S (see e.g. 
chapter 8 in lRiottdl2002n and the CMB temperature anisotropy in 
the Sachs-Wolfe limit AT/T = -0/3 ( .Sachs & Wolfe,, 1967,) . The 
bispectrum of O is defined by 



<0(ki)(l)(k2)(D(k3)) = (2;r)^5o(k, + kz + k3)B^(k„k2,k,). 



(8) 



We will be particularly interested in three bispectral shapes, which 
cover the expec ted shape from a wide range of inflationary models 
( lKomatsull2010l) . They are defined as: 

(i) Local siiape. This is defined by 

B'r'iki , k2,h) = 2A^fl,r' [k'r'kr' + (2 perm)) , (9) 

where the amplitude A is defined by Pt,(k) = Ak'-'^K It may arise 
through a simple Taylor expansion about the Gaussian (linearised) 
perturbation <I)(x) = OgCx) -i- /nlOq, although this is not the most 
general ansatz for <b which gives rise to the local bispectrum. The 
local shape typically arises from super-horizon evolution of the 
curvature perturbation (k <s: a H), which occurs for example in 
some multifield inflation models (Byrnes et al. 2008: EUist on et al.l 
120111) . during modulated reheat ing (Pvali et al. 2004) . in the cur- 
vaton scen ario JLvth et alj|2003l) (the last three models are closely 
connected (^Alabidi et al. 201Cf)), as we ll as the ekpyrotic scena rio 
(^Lehners 2010) and non-local inflation dBarnabv & Clinel2008l) . 
(ii) The equilateral sliape is given by 

B2"\k„k2,h) = 6A\fZ"[-2{hk2k,f"'-'^" 

- [{kikif'^^^ + (2 perm)] 

+ [k["'-'^"kf"'-'^"k^'-' + (5 perm)]) . (10) 

This shape typically arises in models with non-canonical kinetic 
term s, the most studied example being D irac-Bom-Infeld infla- 
tion JSilverstein & Tonsll2004l : lAlishahiha et al. 2004 )). Also var- 
ious other models can produce this shape ( Arkani-Hamed e t al.l 
200 4::ISeerv & Lidsevl[20o3 : IChen et al.|[2007 : Cheung et al. 200"^: 
JLi et'alj|2008l ). 
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(iii) The orthogonal shape is given by 

B%"'''(h,h_,h) = 6AVNL'"'(-8(fcife/t3)'<"'"'"' 
-'i[{kik2f"'-'^^ + (2 perm)] 



+3 [^i- 



-4)/3,2(«,,-4)/3,«,,-4 



■^2 



yt"'"* + (5 perm)]) ,(11) 



dSenatore et alJuOlOf) and it was constructed in order to be ortiiog- 
onal to both the equilateral shape and to a lesser extent the local 
shape. 

The local model is maximised for squeezed triangles fci « fe - k, 
the equilateral model is maximised for equilateral triangles k^ == 
ki - k^ while the orthogonal model receives contributions from 
a broader range of triangles. Another frequently considered shape 
is the enfolded one, which is maximised for "flattened" isosce- 
les triangles ki ^ ki - k^ll, but this can be written as a linear 
combination of the three shapes above. We note that although the 
above three shapes cover many classes of non-Gaussian models, 
there do exist other shapes which cannot be written as a combi- 
nation of the above three shapes including localised or oscillating 
bisp ectra, which may be caused by a fea ture in the inflatons poten- 
tial JChen et al.l2007l : lArroia et al.l201 ih . pa rticle productio n while 
observable modes are crossing the horizon JBamabviBOlOt) (but a 
burst of par ticle production later in inflation generates local non- 
Gaussianitv. lBattefeld et al.ll201 1 ). or an inflaton poten tial with su- 
perimposed oscillations JChen et al.ll2008l : IChenll201 ll) . 

Because /nl for all three shapes is normalised to an equilat- 
eral triangle, but the signal-to-noise is maximised for difl'erent tri- 
angle shapes depending on the configuration it is not surprising 
that the error bars on the th ree shapes are significantly different 
JFergusson & Shellardll2009n . In line with this expectation, we find 
that weak lensing can constrain the local model most tightly and 
the equilateral model least well. When investigating relative mag- 
nitudes between structure formation bispectra and primordial ones, 
we restrict ourselves to the equilateral case, as all primordial bis- 
pectra have equal values fo r this configuration. An altern ative nor- 
malisation was proposed bv lFergusson & Shellara 1 120091) . 

A mild scale dependence of /nl is natural, for bo th the local 
model ( iBvmes et al]|2010l . iBvrnese t al. 2010; Hua n3l2010l), and 
the equilateral and orthogonal models tChenii2005i; JBartolo et al.l 
I2OIOI: iNoUer & Magueii dfl201 ih ISurrage et al .11201 ih . and it may 
be much stronger, e.g. a "step-f unction" which is ze ro on large 
scales and large on small scales JRiotto & SlothluOllf) . Although 
we treat /nl as constant in this paper, the motivation for consid- 
ering scale-dependence is important because weak lensing probes 
smaller scales than the CMB and hence the CMB bounds may not 
ap ply here, see Sec. [4] Observ ational probes have been considered 
in JLo Verde et alJ2008. : .Sefusatti et al.l2009l:IShandera et al.l2oTII : 
lBeckeretal.ll201ll) . 

For converting the bispectrum of the potential fluctuations to 
those of the density field we use the Newtonian Poisson equation 
for ea ch occurence of the potential in the bispectrum iMunshi et al.l 
l20TIh . 

AO = ^^<5 ^ 6(k,a) = ^DAa)(XHkfT(k) <S>(k) (12) 

The horizon entry of each mode is governed by the transfer function 
T(k) and it grows oc D+(a) in the linear regime, such that 

B,{ki,k2,k,,a) = Y\[^D^(")0(Hkd-T{k,)\ B^ikiMM)- (13) 



We choose the normalisation factor A to be consistent for each lin- 
early evolving mode of the density field with our definition of ixg. 

3.2 Non-Gaussianities from structure formation 

Nonlinear processes in structure formation break the homogene- 
ity of the growth equation and generate non-Gaussian features in 
the initially close to Gaussian density field. From Eulerian per- 
turbation theory (see Sahni & Colejl 19951: Bemardeau et al.ll2002l : 
lMatsubarall201ll ; IScoccimarro & Couchmanll2001r) . the first order 
contribution to the bispectrum Bs{ki,k2,ki) (for an introduction, 
see Frv 1984a b) of the density field from nonlinear structure for- 
mation is given by: 

Bs{kuk2,ki,a)= J] D\{a) M(k„kj) P{k,)P{kj), (14) 

ij=l,2.3 

where the classical mode coupling function is M{kj, kj) 

X = kjkjlikjkj) denotes the cosine between the wave vectors A, and 
kj. Due to the fact that P(k,a) grows oc D\{a) in linear structure 
formation, the bispectrum scales with D\{a) in lowest order per- 
turbation theory. In terms of non-Gaussianity parameters and con- 
figuration dependences, structure formation non-Gaussianities are 
strongest for the squeezed configuration because the mode coupling 
function M(A,, kj) assumes the largest values for parallel wave vec- 
tors (with the cosine being one, x = 1), and therefore resembles 
non-Gaussianities of the local type. Their strength in the weak shear 
bispectrum corresponds to an /nl -parameter of (9(10*), i.e. two 
orders of magnitude larger than the primordial non-Gaussianities 
weak lensing can probe. 



4 WEAK GRAVITATIONAL LENSING 

4.1 Convergence spectrum 

The weak lensing convergence k follows from a line-of-sight inte- 
gration weighted with the lensing efiiciency ty,t(Y) (for reviews, see 
iBartelmann & Schneiderll200ll : lBartelmanrJ|2010l) . 



nXH 

K= AxWAx)S 

Jo 



(16) 



and reflects, because of its linearity, all statistical properties of the 
density field d. The weak lensing efficiency is given by 



W.(x)=^-^-^G(j^)x, 



(17) 



with the weighted distance distribution G(x) of the lensed galaxies, 
r"dyq(z)^^^. (18) 



d^' X' 



The spectrum C,t(i) then results from applying Limber's equation 
( lLimbej|l954t) . 



r"" dv 

CM= -^Wl{x)P{k=eix,a). 

Jo X 



(19) 



For the galaxy redshift distribution q(z) we assume a standard 
shape, 

.(.)^.o(|fexp(-(£)'')d. wtth l.|r(|), (20) 
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with the median redshift set to 0.9, as projected for Euclid. 



4.2 Convergence bispectrum 



Similarly as in the case of the weak shear spectrum C^it) we use 
the Limber-equation in the flat-sky approximation, 

fXH d^ 
7 



BAli^ti^li) ■ 



Jo 



W\x) B,(kuk2,ki,a), 



(21) 



with k„ 



f„/x, P 



1 , 2 , 3, for projection o f the flat- 



sky convergence bispectrum 6^ (Schneider et al. 199b; Hul l 19991 : 
iTakada & Jain 2003a b, 2004; Dodelson & Zhang 200^. The 
spherical bispectrum fi^C ^i ,ti,t?,) is related to the flat-sky bispec- 
trum B,{li , ^2, ^3) by CMiralda-Escuddl 1 99 ll ; iKaiseJl 19921) 



BAti^ti^t^) ^ 











An 



B,(t,,t2,ti),{22) 



where 

/A 








^: 



dxPt,{x)Pi^{x)Pt,{x), 



(23) 



denotes the Wigner-37 symbol, which results from integrating over 
three Legendre polynomials P((x) {x = cos 6). The Wigner-3_/ sym- 
bol nulls configurations which would violate the triang le inequality, 
\{i - tj\ < 4 < Yi + lj\ JAbramowitz & Steguj[T97a) . The facto- 
rials in the Wigner-3y symbol are evaluated using the Stirling ap- 
proximation for the F-function, 

r(n + 1) = n! with T{x) == V2Jr exp(-jc) jc^-J , (24) 

for X » 1 JAbramowitz & Stegunlll972h . At this point, it is ap- 
propriate to recall two important issues related to the weak shear 
bispectrum as a line of sight-integrated quantity: The line-of-sight 
integration causes the non-Gaussianities in the convergence to be 
weaker than that of the source field, as a consequence of the cen- 
tral limit theorem, because many uncorrelated lensing effects (if the 
Bo m approximation is invoked and lens-lens coupling is neglected. 



Coorav & Hull2002l;IShapiro & Cooravll2006l: iKrause & Hiratal 



l201(]|;ISchaefer et al.l201 ll) add up to the signal jjeong et al.l201ll) . 
Secondly, the evaluation of the wave vector k = IJx in the source 
field bispectrum B^ generates a mixing of scales when the distance 
X runs over the integration range such that the observed weak lens- 
ing bispectrum is a superposition of density field bispectra of vary- 
ing scale and fixed projected configuration. 

4.3 Properties of the weak lensing bispectrum 

Fig. lAll gives a 3-dimensional visual impression of the three differ- 
ent bispectra as observed by weak shear. The weak shear bispec- 
tra are given as dimensionless bispectra, by multiplication with the 
prefactor (^1^2 ^3)"*'^ ~ ^'*- There are clear differences in the config- 
uration dependence: Local non-Gaussianities provide large ampli- 
tudes for squeezed configurations, i.e. in the corners of the domain 
admissible by the triangle inequality, orthogonal non-Gaussianities 
are largest for folded configurations and the equilateral bispectra 
assume large values if the three multipole orders are equal. 

Fig. [T] illustrates the contribution AC ^(t) I Ax to the spectrum 
and the contribution AB^({, i, i)/Ax to the equilateral bispectrum as 
a function of comoving distance. For the relevant range of multi- 
poles, modes with wave numbers in the range 0.1 ... 1 Mpc/h are 
being probed by weak shear, which are larger than those wave num- 
bers measurable in the primary CMB bispectrum and emphasises 
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Figure 1. Contributions C dC^(f)/dx to the lensing spectrum (solid line), 
and f^dBiiiC, I, i)/dx to the lensing bispectra for both primordial non- 
Gaussianities (dashed line, /nl = 1) and structure formation non- 
Gaussianities (dash-dotted line), for f = 10 (thin lines), £ = 100 (medium 
lines) and { = 1000 (thick lines). The bispectra are plotted for the equilat- 
eral configuration. 



the necessity of mea suring /nl in its scale d ependence, in a similar 
way as advocated bvl Lo Verde et al.l ( l2008li for number counts. 

The skewness parameters S „(£} are defined as the ratio be- 
tween the squared equilateral bispectrum and the cubed spectrum. 



SAO ^ 



cut) ■ 



(25) 



S k{£) is proportional to /nl (and to cr^ to lowest order) and is inde- 
pendent of £1^. It will become relevant for the signal-to-noise ratio 
S(^) (see eqn.[26]in Sect. l5.lt . These skewness-parameters are de- 
picted in Fig.|2]as a function of median redshift Zmd of the lensing 
survey for both the primordial and the structure formation induced 
weak lensing bispectrum. The skewness parameters increase with 
increasing survey depth and in case of primordial non-Gaussianities 
saturate at redshifts of unity, which is an effect of the time evolution 
of the gravitational potential being mapped out by weak lensing, as 
perturbations in the potentials decay in the dark energy-dominated 
phase and are constant in the matter-dominated phase. Structure 
formation non-Gaussianities decrease slightly for deeper reaching 
surveys, which is caused by the fact that the structure formation 
skewness's are building up during f2,„-domination, as they scale 
oc D\/a^ ~ a in contrast to primordial non-Gaussianities in the po- 
tential, which are constant with their scaling oc (£)+(a)/a)' ~ const. 
The plot suggests that with the redshift range probed by Euclid the 
largest-possible skewness's are being observed. 



5 STATISTICS 

5.1 What signal-to-noise ratio can one expect? 

The cumulative signal-to-noise ratio Z(f) for the weak le nsing bis- 
pectrum fi,(fi, ^2, ^3) up to multipole order E is given by ( JHiJ 19991) 



^^w = z z z 



''l=''min''2=''min ''3 = ''m 



cowiiutij^y 



(26) 
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Figure 2. Skewness parameter 5;i;(/') for both primordial non-Gaussianities 
(dashed line, /nl = 1) and structure formation non-Gaussianities (dash- 
dotted line) as a function of the survey depth, for f = 10 (thin lines), t = 100 
(medium lines) and { = 1000 (thick lines). 



Figure 3. Cumulative signal-to-noise ratio S(f)//NL as a function of maxi- 
mum multipole order {, for local (circles), orthogonal (lozenges) and equi- 
lateral (squares) non-Gaussianities. The figure compares the result from dif- 
ferent Monte-Carlo integration routines. 



The summation is carried out with the condition dj ^ £2 ^ ^3 
JTak ada & Jain 2004), such that the covariance becomes 

cov(fi,4,^3) = ^%^^^C(A)C(^2)C(^3), (27) 

,/sky 

where the function A(^i,^2.^3) counts the multiplicity of triangle 
configurations and is defined as 



A(A,^2,^3) = - 



6, i{ = ti = t-i 
2, ti = ij for i 9^ ;• 
1, t^i^e^tt-i* i\ 



(28) 



/sky denotes the fraction of the observed sky and is set to /s^y = 1/2 
for Euclid. The observed spectra 

al 
C(l) = C(£) + — , (29) 

n 

with the number density of ellipticity measurements per steradian 
n, which is set to 40 galaxies per squared arcminute, corresponding 
to the projected Euclid performance. Instead of a direct summation 
over £\,£2, £3 we use a Monte-Carlo integration technique and con- 
sider the evaluation of eqn. | |26! > as a three-dimensional integration, 
for which we use publicly available CUBA-library (Hahn 2005) . 

The cumulative signal-to-noise ratio !,(£) for a measurement 
of the weak shear bispectrum is depicted in Fig.[3]as a function of t 
and for all three non-Gaussianity types. As the signal strength E is 
proportional to the non-Gaussianity parameter /ml, it is convenient 
to plot the ratio S(f)//NL. The plot suggests that weak lensing bis- 
pectra sourced by primordial non-Gaussianities could only be mea- 
sured with Euclid for /nl significantly larger than 100. Orthogonal 
bispectra are weaker by a factor of 3 compared to local bispectra, 
and equilateral bispectra generate the weakest signal, being a factor 
of 8 weaker than local non-Gaussianities. The four different MC- 
integration algorithms agree well in their results for !,(£}, and when 
the number of sampling points is chosen to be oc £, the algorithms 
retain accuracies, indicating that the adaptive algorithms take ac- 
count of the symmetry properties of the integrand. At f = 1000, is 
is sufficient to compute 0(l(fi) samples, which reduces the num- 
ber of evaluations of the integrand by a factor of lO'' compared to 
the exact evaluation, which would require (9(10') evaluations, for 



precisions on the sub-percent level. We will restrict ^-space inte- 
grations to multipoles i 1000, because the increase in signal when 
extending the f-range is marginal for primordial non-Gaussianities 
and additionally, it helps to avo id scales influenced by bary onic 
physics and intrinsic alignments JSemboloni et alj2008ll201 11) . 

The configuration dependence of the contribution to the in- 
tegrated signal is given in Fig. lAll in the appendix for the three 
bispectrum types considered here. The panels show the weak shear 
bispectrum in units of the noise, B;,(Ci,{2, ^3)/ '^cav{l\,t2,ti) as a 
function of ti, £2 and ^3. While the local non-Gaussianity provides 
the largest contributions for squeezed configurations, the orthogo- 
nal non-Gaussianity shows a much uniformer contribution through- 
out ^-space, and the equilateral non-Gaussianity is only providing 
significant amplitudes for very small values of £. This behavior is 
reflected in in the cumulative signal-to-noise ratio, as shown in 
Fig. [3] From the signal-to-noise ratios one can already estimate 
the accuracy for a measurement of /nl- The conditional Cramer- 
Rao bounds (T/^j^ = 1/ Y^aJlaJl "'i ^^^ non-Gaussianity parameter 
(which are at the same time the non-Gaussianities required to gen- 
erate a signal of unity) are cryj^^ = 200 (local), CTf^^^ = 575 (orthog- 
onal) and (Tfj^j^ = 1628 (equilateral) which is significantly weaker 
than other probes such as the primary CIVI B, due to the Gaus sianis- 
ing effect of the line-of-sight integrations jjeong et al.ll201lh . 



5.2 Would one misestimate /nl using the wrong bispectrum? 

The ^^-functional constructed for measuring the noise-weighted 
mismatch between the true bispectrum B[ and the wrongly assumed 
bispectrum 6JJ' for interpreting the data reads: 



X' 



z z z 



[aB:(£i,£2J3)-Bi(£iJ2J3)] 

COV(A,^2,^3) 



(30) 



and yields the best fitting a from the minimisation dx^ Ida = 0. 

The variable a measures the ratio between the wrongly in- 
ferred non-Gaussianity parameter /nl and the true value, and is 
given by Fig. |4] as a function of maximum multipole order con- 
sidered in the integration of the ^^-functional. For weak signals 
this ratio is very close to unity, and differences emerge when the 



6 B.M. Schdfer, A. Grassi, M. Gerstenlauer, C.T. Byrnes 
















' 

















- local 

























- ortho 




















o 


n 


equil 
local - 
ortho - 
equil ■ 


0"^ 










^*^ 




^ 
















^ 




^ 




















\ \ 




\ 








0"^ 












1 \ fa 


\ 
\ 


\ 








0"^ 












\ \ 




\ 
\ 
\ 


\ 






■* 
















1 
1 




|l 






1 i '' 



multipole order £ 



10 10 10 10 

non-Gaussianity parameter /nl 



Figure 4. The ratio a{t) of inferred /NL-value to tlie true /nl -value as a 
function of maximum multipole order t. Tlie true non-Gaussianity model 
is indicated by the line style, whereas the wrongly chosen non-Gaussianity 
model is given by the marker style: local (circles, solid lines), orthogonal 
(lozenges, dashed lines) and equilateral (squares, dash-dotted lines) non- 
Gaussianities. 



Figure 5. Probability <?(/nl) of obtaining data as extreme as a wrong fit 
to the weak lensing bispectrum by choosing the wrong non-Gaussianity 
model (indicated by the line style) to data (indicated by the marker style): 
local (circles, solid lines), orthogonal (lozenges, dashed lines) and equi- 
lateral (squares, dash-dotted lines) non-Gaussianities. The horizontal lines 
indicate 1, 2, 3, Act confidence intervals and C has been set to 1000. 



integration is carried out to larger multipoles, and the signal be- 
comes stronger. Misestimations in /nl up to half an order of mag- 
nitude appear possible, including wrong signs for /NL-estimates. 
Most combinations of fij^ and B" yield very small values for the 
estimated /NL-parameter (equivalently, a ^ —\) when choosing the 
wrong non-Gaussianity. Again, the evaluations necessary for de- 
termining a are carried out as an MC-integration with the CUBA- 
library ( lHahnll2005h . 

5.3 Would one notice fitting tlie wrong bispectrum? 

Now, the question appears if one would notice the assumption of a 
wrong primordial bispectrum when fitting for the non-Gaussianity 
parameter /nl. This can be quantified by the probability q of ob- 
taining data more extreme than the one at hand. This probability q 
(Fisher's p-value) is given as a function of the true /nl under the 
assumption of a Gaussian likelihood, which is well justified given 
the very large number of degrees of freedom (although doubts have 
been raised about how accurate this is, see ISmithet al.ll201ll) . As 
shown in Fig.|5]for all combinations between true non-Gaussianity 
types and wrongly fitted non-Gaussianity models, this probability 
drops very rapidly towards very small numbers for /NL-values of a 
few hundred, indicating that it would be very difficult to reconcile 
non-Gaussianities of that strength with observations if the wrong 
non-Gaussianity model had been chosen. For /NL-values smaller 
than 100 the signal is so weak that no significant discrepancies be- 
tween data and model appear, for any type of non-Gaussianity. 
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Figure 6. lo", 2o" and Sir-constraints on the non-Gaussianity parameter /nl 
and the matter density n,„, for the local non-Gaussianity (solid lines) and 
the orthogonal non-Gaussianity (dotted lines) at f = 1000, from a numerical 
MC-integration of the Fisher-matrix. 



5.4 Do parameter constraints depend on non-Gaussianity? 

The Fisher-matrix formalism JTegmark et alJI 19971) is widely used 
in cosmology for deriving parameter forecasts, and requires in the 
case of the bispectrum as the signal-to-noise ratio the summation 
over all triangle configurations: 






dB, 



1 



as. 



dx,, co\{ti,l2,li] dxy 



(31) 



The explicit summation can be replaced by a d'' ^-integration, which 
can be carried out using the MC-technique outlined in Sect. 15.11 
Resulting simultaneous constraints on fi„, and /nl from the weak 
shear bispectrum sourced only by primordial non-Gaussianities 
with no other priors are given in Fig. |6] for local and orthogonal 
models at a reference /nl = 1000. The equilateral bispectrum does 
not constrain the parameter pair in a meaningful way due to the 
weak signal. 
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6 SYSTEMATICS DUE TO STRUCTURE FORMATION 

6.1 Can one subtract the structure formation bispectrum? 

Naturally, the small primordial non-Gaussianities are superseded 
by much stronger non-Gaussianities due to nonlinearities in the 
cosmic structure formation processes, which affects the measur- 
ability of /nl- There exists an accurate description of the struc- 
ture formation bispectrum provided by Eulerian p erturbation theory 
on th e scales of interest (compare Sect. 13.21 and lSemardeau et al.l 
l2002h . if the cosmology is known - but there are always uncer- 
tainties in the cosmological parameter set, which would result in 
an uncertainty in predicting the structure formation bispectrum. If 
the structure formation bispectrum is not properly subtracted from 
the observed bispectrum, there will be errors in the estimation of 
the non-Gaussianity parameter /nl for the primordial bispectrum, 
which can be quantified with the;j'^-functional. 



X 



I. I. Z 

^l~^inin ^2— '■min 0~'mi 



COV(A,^2,^3) 



(32) 



describing the fit of a primordial bispectrum Bj^ to data AB^, 

AB,(A , ^2, ^3) = B[(ft J2Ji) + Sf (A , i2Ji)-fff\l, , ^2, ^3)(33) 

which contain the true primordial bispectrum B'^ itself, the very 
large structure formation bispectrum B'-^^ for the true cosmology, 
from which the structure formation bispectrum B"-^^ has been sub- 
tracted, possibly incompletely, by assuming the wrong cosmology. 
Derivation dx' Ida = yields the best fitting a, which is related 
to the misestimated f^^^ = a/NL and the deviation from the true 
non-Gaussianity S = a - I. 

Distributions p(S)d6 have been derived for all bispectrum 
types by drawing 10' samples from a Gaussian likelihood for the 
parameters SI,,,, erg, h, iis and w of a standard spatially flat dark 
energy model. The covariance matrix has been constructed using 
the icosmo resource for the Euclid weak lensing and BAG data 
JRefr egier et al. 201 1), and provides an excellent prior on the cos- 
mological parameters. By this sampling process of S it is possible 
to propagate the entire uncertainty in the cosmological parameter 
set onto the estimate of /nl. As shown by Fig.l?] the resulting dis- 
tribution is very close to Gaussian, with zero mean and standard 
deviations of ixy^^ = 119 (local), cTfj^^ = 372 (orthogonal) and 
''"/nl =511 (equilateral), which is similar to the statistical uncer- 
tainty of measuring /ml and thus constitutes a serious error. The 
width of the distributions are independent of the true value /nl, 
and the relative error (5//nl scales oc 1//nl. Misestimates of that 
magnitude make it very difficult to assign a primordial origin to a 
non-zero residual bispectrum, given the current bounds on /nl. All 
integrations were computed up to ^ = 1000. 
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Figure 7. Distiibutions p(d)<i6 of the bias S between the inferred non- 
Gaussianity parameter and the true parameter /nl if the structure formation 
bispectrum is not completely removed, due to uncertainties in the cosmo- 
logical model, for local (circles, solid line), orthogonal (lozenges, dashed 
line) and equilateral (squares, dash-dotted line) non-Gaussianities. As pri- 
ors, Euclid weak lensing and baryon acoustic oscillations were used. 
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Figure 8. Distributions p(S)dS of the bias 6 with an enhanced prior for the 
cosmological model with constraints from weak lensing, baiyon acoustic 
oscillations (both Euclid) and CMB temperature and polarisation spectra 
(Planck), for local (circles, solid line), orthogonal (lozenges, dashed line) 
and equilateral (squares, dash-dotted line) non-Gaussianities. 



6.2 What happens if a better prior is available? 

The uncertainty in predicting the weak shear bispectrum generated 
by nonlinear structure formation can be reduced if a stronger prior 
on the cosmological parameter set is available or if the complexity 
of the model is reduced, e.g. if the wCDM dark energy cosmol- 
ogy would be replaced by the simpler ACDM cosmology. Fig. [8] 
illustrates the distributions p(6)dS of the difference S between the 
inferred non-Gaussianity parameter and the true parameter, if the 
contamination of the bispectrum is computed by drawing 10' sam- 
ple wCDM cosmologies from a Gaussian likelihood whose covari- 
ance matrix incorporates constraints from Euclid weak shear spec- 
tra, Euclid baryon acoustic oscillations and in addition Planck's 



constraints from the observation of primary CMB temperature and 
polarisation spectra. In comparison to the distributions shown in 
Fig. |7] the width is now much reduced, by about a factor of 4, 
allowing measurements down to smaller values for /nl. The spe- 
cific uncertainties are CTf^^ = 29 (local), (Tj^j^ = 98 (orthogo- 
nal) and (Tf^^ = 149 (equilateral), all stated as standard devia- 
tions of the distribution p(d)d6. These uncertainties are below cur- 
rent bounds on /nl and are small enough for studies of primordial 
non-Gaussianities. Again, all integrations were carried out up to 

e = 1000. 
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7 SUMMARY 

The topic of this paper are measurements of primordial bispec- 
tra in weak shear data from Euclid, comparisons between differ- 
ent types of non-Gaussianity configuration dependences, statisti- 
cal questions concerning the inference of the non-Gaussianity pa- 
rameter /nl and the removal of the much stronger structure forma- 
tion induced bispectrum. Although not as sensitive as observations 
of the CMB-bispectrum or the galaxy bispectrum for scale-free 
non-Gaussianities, weak lensing can place useful independent con- 
straints on non-Gaussianities, in particular on smaller scales where 
CMB bounds might not apply. It is less prone to systematics than 
other large-scale structure probes and provides a direct linear map- 
ping of the density field, which conserves its statistical properties. 

(i) Primordial non-Gaussianities provide a rather weak signal 
in the weak shear bispectrum (because of the Gaussianisin g ef- 
fect of the line-of-sight integration, see Ueong et alj uOllh . and 
signal-to-noise ratios of order unity can only be expected for 
/nl = 200,575,1628 for local, orthogonal and equilateral non- 
Gaussianities, respectively, where this measurement is most sen- 
sitive to scales 0.1 ... 1 (Mpc//i)"'. These bounds are weaker than 
those from e.g. observations of the CMB bispectrum, but will serve 
nevertheless for cross validation, in particular given the absence 
of strong systematics is weak shear data, or as bounds on scale- 
dependent non-Gaussianity, because weak lensing maps out scales 
which are not constrained by the CMB and is sensitive to scales 
probed by number counts. 

(ii) Configuration space integrations can be very efficiently car- 
ried out by Monte-Carlo integration schemes, at a fraction of the 
computational cost. Computations of the signal-to-noise ratio, of 
^--functionals or of the Fisher-matrix F,„, can be done with accu- 
racies below a percent with (9(10^) evaluations instead of (9(10') 
evaluations for the direct sum over li , ti and tj,. V ery good results 
were obtained with the CUBA library JHahnll2005h . 

(iii) Fitting the wrong bispectrum type to data yields serious 
misestimates in the non-Gaussianity parameter /nl- Depending on 
the combination of true and false model there are two cases: either 
the estimated /nl becomes very small, or the estimate for /nl is 
a factor of ~ ±3 too large. When looking at numerical values for 
the ^^ -functional, one would notice strong discrepancies between 
data and model when fitting the wrong non-Gaussianity type from 
values of /nl of a few hundred on. 

(iv) The much stronger structure formation bispectrum can be 
subtracted with a prediction of its bispectrum from perturbation 
theory if the cosmology is known precisely enough. Propagating 
the uncertainty in the cosmological parameter set onto the misesti- 
mation of /nl if the structure formation bispectrum is not correctly 
subtracted yielded typical uncertainties of 29, 98 and 149 for lo- 
cal, orthogonal and equilateral non-Gaussianities, much less than 
the statistical accuracy. As a prior on the cosmological parameters 
we assumed a Gaussian likelihood for a wCDM model combining 
Euclid's weak shear with baryon acoustic oscillations and Planck's 
observations of primary CMB anisotropics. Similar ratios between 
the numerical valu e of /nl and the s tandard cosmological parame- 
ters were found bv lPace et al.l <201lh . 

Many of our investigations can be straightforwardly generalised 
to other probes of large-scale structure statistics. We intend to 
generalise our investigations to highe r poly spectra and t o apply 
ideas from Bayesian model selection ( lTrottall2007l 120080 for as- 
signing probabilities to the problem of choosing the correct non- 
Gaussianity type. 
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APPENDIX A: CONFIGURATION DEPENDENCE 

Fig. lAll compares the configuration dependence of the bispectrum 
and of the signal strength in a weak lensing experiment. As a 
representation, we chose to plot the dimensionless weak conver- 
gence bispectrum (^i^2^3)'''^S,(^ii^2.^3) and the convergence bis- 
pectrum in units of the noise, B,(€i, ^2, ^3)/ '^cov(£i,£2,fi), which 
when added in quadrature yields the signal-to-noise ratio. The fac- 
tor (fi^2^3)'''^ ~ ^* makes the angular bispectra dimensionless. 

This paper has been typeset from a TgX/ ffl^ file prepared by the 
author. 
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Figure Al. Configuration dependence (fi^2^3) ^it(^l.^2,^3) (first column) and signal-to-noise ratio B/j(^i , {2, ^3)/ iJcov({i , {2, ^i) (second column) of the 
weak lensing bispectrum, for local (first row), oilhogonal (second row) and equilateral (third row) non-Gaussianities. The size of the blobs and their colour is 
proportional to the bispectrum, where a correct relative normalisation in the columns is maintained. Configurations outside the grey bounding planes violate 
the triangle inequality. 



